Quantum Dynamics of a Nanomagnet driven by Spin-Polarized Current 
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A quantum theory of magnetization dynamics of a nanomagnet as a sequence of scatterings of each 
electron spin with the macrospin state of the magnetization results in each encounter a probability 
distribution of the magnetization recoil state associated with each outgoing state of the electron. 
The quantum trajectory of the magnetization contains the average motion tending in the large spin 
limit to the semi-classical results of spin transfer torque and the fluctuations giving rise to a quantum 
magnetization noise and an additional noise traceable to the current noise. 
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Introduction The spin transfer torque (STT) [l|, Q is 
both a fundamental process in magnetization dynamics 
and important to spintronics in memory and informa- 
tion processing integration. The magnetization dynam- 
ics driven by spin-polarized current via STT has made 
great progress by experiments and by a "semiclassical" 
theory of treating the magnetization dynamics classically 
and the polarized current transport and the spin waves 
quantum mechanically (see a series of reviews introduced 
by Ref. d). 

Recent development in fast time-resolved measure- 
ments and coherent control have made possi- 
ble studies of magnetization dynamics and fluctuations 
not masked by the inhomogeneity effects of the measure- 
ment and the prospect of precision control. These recent 
experiments including Ref. [8[, by showing the stochas- 
tic nature of the magnetization dynamics at short time 
intervals, highlight three important aspects of the sub- 
ject. First is the stochastic motion which may lead to 
an understanding of noise, a fundamental issue in mag- 
netization dynamics 9]. Second, in interconnected sys- 
tems at the nanoscale, the particulate nature of the cur- 
rent electrons imprints shot noise on the magnetization 
through the quantum scattering process USES. Third, 
the thermal fluctuations are important because, for ex- 
ample, the thermal noise power is comparable to the mi- 
crowave power generated by the magnetization preces- 
sion [12]. The measurements seem to be satisfactorily 
treated by the semiclassical theory including micromag- 
netics. However, the common treatment of angular mo- 
mentum transfer by the magnet as an immovable scat- 
tering potential is shown below by order of magnitude 
estimate to be inadequate for accurate computation of 
mean displacement and fluctuations in a nanomagnet. 
The origin of the stochastic motion arises from the as- 
sumption of randomness in either the current or the mag- 
netization. The spin waves seem grafted on rather arising 
naturally out of the vibration of the spins which consti- 
tute the moving magnetization. These points motivate 
us to question if inclusion of the quantum dynamics of 
the magnetization may not only provide more accurate 
connection between the current electrons and the magne- 



tization but also provide qualitative understanding of the 
stochastic dynamics as well as bringing out the prospect 
of coherent motion for precision control. By a quantum 
theory of the nanomagnet which gives its mean dynamics 
and fluctuations through scattering between the current 
electrons and the movable magnet, we hope to illustrate 
the fundamental aspect of our approach to STT. The de- 
velopment of quantum optics after the laser operation 
was understood by semiclassical theory provides perhaps 
an optimistic historical guide for the development of co- 
herent magnetization dynamics after the successes of the 
semiclassical STT theory. 

Conditional Scattering and Stochastic Schrddinger 
Equation We take the magnetization of a single-domain 
nanomagnet to be represented by a coherent state of 
the local spins [13[ and view its dynamics in the spin- 
polarized current as a sequence of scattering by individ- 
ual electrons. Since each scattering produces a number 
of outgoing macro-spin states, the serial scatterings cre- 
ate a Monte Carlo trajectory of the macro-spin states 
which is governed by a stochastic Schodinger equation 
14| . Repeated solutions form an ensemble of quantum 



trajectories for the magnetization from which the mean 
trajectory and the fluctuations can be obtained. Since 
the outcomes of the macro-spin states are paired with 
the scattered spin states of the current electron to form 
an entangled state, the scattered electron state may be 
viewed as conditioned on the corresponding macro-spin 
state. After the first scattering, the electron recedes and 
the magnetization is left in a reduced density matrix of 
the outgoing macro-spin states. The second electron from 
the current may be viewed as the recipient of one of 
the macro-spin states. This step is equivalent to mea- 
surement or decoherence 15], and the magnetization un- 
dergoes a Brownian motion and not a quantum random 
walk. The magnetization noise then has three sources, 
an intrinsic one due to the uncertainty of the quantum 
state of the magnetization, two extrinsic ones due to the 
entangled state after each scattering and due to the par- 
ticulate nature of each colliding electron with the coher- 
ent macro-spin state. The last forms a channel for the 
mutual effects of the current noise and the magnetization 



noise. 

We illustrate the quantum effects with a simple scat- 
tering model of the nanomagnet by a regular sequence 
of electrons and then layer the complexities of the real 
current 16]. We study here the dynamics of the magneti- 
zation as an angular momentum coherent state |J, 0,3?) 
(l3j . The effect of the spin waves as small deviations from 
the coherent state is being investigated. J is the quantum 
number of the total spin J of the localized d-electrons in 
the nanomagnet and is of the order of 10 6 for typical 
nanomagnets. The angles and <f> characterize the di- 
rection of the total spin of the nanomagnet. The current 
electron has the incoming state |fc,t) 5 its wave vector k 
impinging normally on a ferromagnet film in the x = 
plane and its spin along the z-axis. The scattering po- 
tential, <5(x)[AoJo + As • J], is a potential of fixed position 
but includes a spin independent term indicated by the 
identity operator J of the total spin state space and the 
exchange term between the current electron spin and the 
total local spin. The <5-potential models the FM film of 
thickness of several nanometer, and the parameters Ao 
and A are the interaction strength between the current 
electron and one ci-electron, estimated from two simple 
exchange-split potentials. 

After a scattering event by one current electron, the 
total state of the system of the nanomagnet and the cur- 
rent electron is entangled In the Berger Q basis 
states {| ± k, ±)} for the current electron, in which the z 
axis is rotated in the z-3 plane to the magnetization di- 
rection before the scattering, the conditional macro-spin 
states associated with the outgoing single electron states 
are coherent states with an error of the order 1/J 2 , 
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\V)=J29i\J,®u®i)\ki,*i), (1) 

where Span(|A;i, <n)) = (| - fc,+), \k,+), | - k, -), \k, -)). 
The Berger basis is therefore preferable to the original 
spin basis of the current electron of spin up and down 
along the z -axis, and only the incoming electron state 
is shown in the original basis. The Berger basis may 
be viewed as constituting a measurement basis for the 
current electron with the outcomes associated with the 
probabilities Gi — \gi\ 2 ,i — 1 — 4. A conditional state 
of the macro-spin state of the magnetization is deter- 
mined by the scattered current electron state, as shown 
by the Feynman diagrams in Fig.[T] Note that, after scat- 
tering, the conditional macro-spin states associated with 
the electron states | ± k, — ) suffer only a phase change be- 
cause the flip-flop terms s+J_ + sLJ + in the Hcisenberg 
exchange interaction do not connect the state \k, +)\J, J) 
to \k, —)\J,m) for any m < J. 

The evolution of the quantum state of magnetization 
from |J, 0,<f>) to |J, 0i,$j) stochastically by one spin- 
polarized current electron may be represented by a ro- 
tation, | J, 0',$') = e- ifl ' Jl? |J,0,$), with the rotation 
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|-k,+> |J,0r,<J> ,■) |k,+> |J,0„<JO |-k,-> |J,0,<t>> |k,-> |J,0,O> 




k,t> |J,0,<t>> |k,t> |J,0,O> |k,'> |J,0,O> |k,t> |J,0,<t>> 

FIG. 1. Four possible states of the nanomagnet in the Berger 
basis {| ± k, ±)} after scattering. The straight lines repre- 
sent the spin-polarized electron and the wavy lines represent 
the nanomagnet states. The macro-spin state J, Oi, <3>i) with 
suffices i =1 — 4 define the association with the electron spin 
states in the Berger basis. Gi are the conditional probabili- 
ties. 

axis n = (— sin ip, cos ip, 0) and the rotation angle d as 
the random variables with probability outcomes deter- 
mined by Eq. (flj. Since the angle d is of the order of 
1/J and the time interval r between two successive elec- 
trons is small compared with the time scale of magnetiza- 
tion dynamics, the stochastic rotation of the nanomagnet 
may be treated as a continuous rotation governed by the 
stochastic Schrodinger equation, 

z^|J,0,$) = W n-J|J,0,<f), (2) 

under the effective magnetic field due to the polarized 
current electrons, producing precession frequency, oj = 
d/r about the axis n. 

Dynamics and Fluctuation of Magnetization, Electric 
Current and Current Noise To study the switching be- 
havior of nanomagnet at zero magnetic field, we com- 
puted an ensemble of 500 runs of Monte Carlo steps, each 
with sequential scatterings by N e = 1.5 x 10 7 electrons 
with the macro-spin J = 10 6 . To allow for the three di- 
mension nature of the electron wavevector distribution in 
a Fermi sphere of radius kp in the normal metal model, 
the normal component to the interface k obeys the dis- 
tribution function f(k) = 2k/k F for k G [0,fc F ][I3|. In 
our simulation, the wavevector k of each injected elec- 
tron is generated randomly according to f{k). The scat- 
tered state in each step is randomly selected according to 
the probabilities Gi , associated with the macro-spin state 
| J, 0j, <F;) and output electron state \ki, ai). The effect of 
STT comes from the random encounter by each incom- 
ing electron of a choice of the macro-spin states left by 
the previous electron, giving rise of a random walk of the 
magnetization. Fig. [2] shows the direction of the magne- 
tization J averaged over the ensemble of trajectories and 
its uncertainties in the Cartesian components AJ Q . 500 
runs were found to be sufficient to stabilize the ratio of 
the uncertainties to the average. 

Fig. 12a) shows that the switching motion also con- 
tains a precession about the spin polarization direction 
of the injected electrons, which reproduced the semiclas- 
sical results of being driven by M x S and M x (M x S) 
for the magnetization M and current spin polarization S. 
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FIG. 2. (color online). Time evolution of magnetization and 
its noise, (a) The average trajectory of the macro-spin state 
under STT and its projections onto the Cartesian planes. 
The Cartesian components J a ,a = x,y,z are normalized in 
units of J. (b) The time evolution of the uncertainties of the 
macro-spin components and its projections onto the Cartesian 
planes. AJ a are normalized in units of J. The current elec- 
tron number is N e — 1.5 x 10 7 and the macro-spin J = 10 6 . 
The initial orientation is (0o,"3?o) = (3,0.5) rad. The elec- 
tron wave vector obey the probability function f(k) = 2k/k F 
for k £ [0, fcf] with kp = 13.6 nm -1 , and the spin polariza- 
tion vector is taken S = (0,0,1). Ao and A are determined 
by the spin-dependent potential A+ = 1.3 V, A_ = 0.1 V, 
and layer thickness d — 3 nm. The switching time between 
points A and C is estimated as 1.3 ns for bias V = 1 mV and 
cross-section area A = 10 4 nm 2 , by Eq. ©. 



It is also found that the number of electrons needed for 
switching is of the order of J. For typical experiments 
17] . the electric current is of the order of mA and the 
switching time is of the order of ns, both consistent with 
the electron number N e we used. Since the input pa- 
rameters to our model, namely the Heisenberg exchange 
energy and the current, are of the order of the physical 
properties of real systems, we see the ball-park agree- 
ment of the switching time as an encouraging sign for the 
quantum approach so that the noise effects from the the- 
ory may be worth being included with the thermal noise 
in interpretation of experimental measurements. With- 
out the thermal noise, the fluctuations AJ Q , shown in 
Fig. [2jb) in time order from point A to C, first increase 
till the midpoint B and then decrease to C. The ini- 
tial state of nanomagnet is a pure angular momentum 
coherent state | J, 0, $) with only intrinsic magnetization 
noise (AJ a ~ \/~J). Scattering by current electrons leads 
to stochastic motion which accounts for the rise of the 
extrinsic magnetization noise. As all the trajectories con- 
verge to the state | J, 0, 0), the magnetization noise finally 
decreases to the intrinsic quantum fluctuation. 

The injected current electron noise can come from ei- 
ther its wave vector distribution f(k) or the uncertainty 
of its spin state. To identify their effects on the magne- 
tization dynamics of nanomagnet, we include them sep- 
arately in the simulations and compared the case with 
no injected electric current noise. As shown in Fig. 02 
the inhomogeneous effect of k has less effect on the time- 
evolution of average value (0) and (<£>) and their fluc- 
tuations A0 and A$; while the electron spin noise can 
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FIG. 3. (color online). Effects of injected current electron 
noise on the magnetization dynamics of nanomagnet. (a) 
Time evolution of the average value of 0. (b) Time evolution 
of the average value of $. (c) Time evolution of the fluctua- 
tion of ©. (d) Time evolution of the fluctuation of $. Three 
cases are considered here: the fully polarized electrons with a 
single wave vector 2k f/ '3 (solid line); the fully polarized elec- 
trons obey wave vector distribution function f(k) — 2k/ k F 
for k £ [0, fop] (dash line); partially polarized electrons with 
a single wave vector 2k F /S (dotted line). The spin polar- 
ized vector is S = (0,0, P). The current electron number is 
N e = 3.0 x 10 7 , and the other simulation parameters are the 
same as in Fig. [2] 



change the magnetization dynamics significantly. The 
magnetization switch time for spin polarization vector 
S = (0,0,0.5) is about twice of that for the fully- 
polarized electrons. The magnetization noise is also ob- 
viously enhanced, although its main contribution is still 
from the scattering process shown in Fig. [TJ 

The electric current I and current shot noise 5 in a 
nanomagnetic junction are given as [18[ 

9 N c 9 N c 



I =^n V E T n,S = 2e— V5>(l-T w)l (3) 

n— 1 n— 1 



where T n is the transmission probability of the electron 
with mixed spin states reduced from the entangled state 
of Eq. ([IJ. iV c is the number of transport channels, esti- 
mated as N c ~ Akp. The channel sums can be replaced, 
respectively, by the ensemble average NcS^ and variance 
N c (A6u) 2 of the current electrons, where 8^ is a random 
variable given the value 1 or for a transmitted or re- 
flected electron. By Eq. ((3|), the time interval r between 
two successive injected electrons is given in terms of the 
current as r = 2nh/N c eV . Then the number of scatter- 
ing events in the simulations can be converted to the time 
scale of magnetization switch. Fig.|4]shows the calculated 
electron current / and current noise S associated with the 
magnetization dynamics in Fig. [3] The electric current 
decreases during the magnetization switch, because the 
electrons are scattered by spin-dependent potential bar- 
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FIG. 4. (color online). Output electric current and current 
noise, (a) Current during the switching of magnetization in 
Fig- El (b) Current shot noise during the switching of mag- 
netization in Fig. [3] 



riers of the nanomagnet, A_ = 0.1 V and A + = 1.3 V, at 
the starting point and final point respectively. The calcu- 
lated current shot noise S rise and drop in sequence, with 
the maximum amplitude corresponding to the transmis- 
sion probability T — 0.5 for electrons. The calculated 
current amplitude and switching time are comparable to 
the experiment results [l7|. The amplitude of shot noise 
of the order of nA 2 /Hz is also consistent with the noise 
measurement |19j . We found that the change of current 
and current noise during the magnetization switch de- 
crease with the reduction of current electron spin polar- 
ization. 

The quantum noise generated in magnetization and 
electric current due to the scattering can be veiled by 
the thermal noise in the system at high temperature. 
The current shot noise will become important when 
eV > kT[HI Eli where K is the Boltzmann constant 
and T is the temperature. For the bias V = 1 mV, 
this gives T < 11 K, and the shot noise was ob- 
served experimentally at 12 K[19j|. The magnetization 
thermal noise can be characterized by the correlator 
(h(t)h(t')) = (2a^hMV)6(t - t') for the random ef- 
fective field h(t)^, [ill, Hi, where a is the Gilbert 
damping coefficient, 7 is the gyromagnetic ratio, M is 
the magnetization, and V is the volume for the nano- 
magnet. Then the thermal fluctuation for the rotation 
angle = 7 f* +T h(t')dt' in the time interval r will 
be {5{) 2 )therm = 2oiqkTt / JH, here we used the relation 
~fjh = MV. On the other hand, the quantum noise 
based on Eq. © gives {8-d 2 ) quant ~ 0{l/J 2 ). Thus the 
characteristic temperature that the quantum noise will 
be comparable to the thermal noise can be roughly esti- 
mated by the relation 



(4) 



smaller Gilbert damping coefficient, smaller nanomagnet, 
or shorter magnetization switch time. 

Conclusion The quantum trajectory method we have 
developed not only reproduced the semiclassical results 
of STT, but also revealed the generation of quantum 
noise and the noise transfer between the nanomagnet and 
current. More quantum phenomena in the STT-related 
physics are anticipated in further studies based on this 
quantum picture. 
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